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Several properties of the solar interior are determined with a very high 
accuracy, which in some cases is comparable to that achieved in the deter- 
mination of the Newton constant Gm- We find that the present uncertainty 
AGn/Gn = ±1.5 • 10 -3 has significant effects on the profile of density and 
■ pressure, however it has negligible influence on the solar properties which can 

be measured by means of helioseismology and 8 B neutrinos. Our result do not 
support recent claims that observational solar data can be used to determine 
^3 ' tire value of Gat with an accuracy of few part in 10 . Present data cannot 

constrain Gjy to much better than 10~ 2 . 

^ I. INTRODUCTION 

o ■ 

In the last ten years, our observational knowledge of the solar interior has progressed 
enormously. By means of helioseismic data it has become possible to derive the sound 
speed in the solar interior with accuracy of about one part per thousand By the 

same method, it has been possible to deduce important properties of the convective enve- 
lope. The photospheric Helium fraction Y p h and depth of the convective envelope Rb have 
been determined with accuracy of about one per cent and one per thousand respectively, 
following the pioneering papers by H and [Q. The measurement of the neutrino flux from 
Boron decay, obtained by combining SNO and Super-Kamiokande data has provided 

a determination of the temperature T c near the solar center with accuracy of about one per 
cent 0. 

It is impressive that all the predictions of the Standard Solar Model (SSM) have been 
confirmed by these accurate tests, see e.g. [0,0 and the first rows of table | for a summary 
of the available information. 

The SSM, like any stellar evolution calculation, depends on several parameters. In this 
respect the Newton constant Gn plays an important role, since stellar evolution results 
from the equilibrium between gravitation and other interactions. The sensitivity of stellar 
evolution to modifications of was first stressed by Teller |K|. By means of a homology 



argument he demonstrated that the stellar luminosity is L & oc G 7 n Mq, so that even a ten 
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per cent variation of Gn from its standard value would imply that life on Earth cannot be 
sustained! 

Among the "fundamental parameters" of nature Gn is the most difficult to 



measure. In fact, it is determined with relatively poor accuracy and the situation of the 
field can be summarized by observing that the CODATA-98 [I^| value, Gn = 6.673(10) 10 -11 
m 3 s _2 Kg -1 includes a relative accuracy AGn/Gn = ±1.5 ■ 10 -3 which is a factor ten larger 
than that estimated in 1986 []TJJ. Recent experiments have claimed accuracy of 

about 10~ 5 however the disagreement between individual results is at the level of 10~ 3 . 

Since several solar properties are now observationally determined with an accuracy of 
order 10 -3 , we shall address the following questions: 

1) What is the uncertainty on SSM predictions induced by the uncertainty o/Gjv? 

2) Can one exploit the available accurate observations of the solar interior, i.e. helio- 
seismic data and Boron neutrino flux, for obtaining a determination of Gn with accuracy 
comparable or better than that of laboratory measurements? 

This last question, which is particularly interesting, was recently raised by Lopes and 



Silk H7I 



We shall present solar models calculated for different value of Gn, deriving the predictions 
to be compared with helioseimic and neutrino data. We shall find that solar sound speeds and 
temperature are actually very weakly dependent onGjy and we shall provide an explanation 
of this apparently puzzling result. 



II. SOLAR MODELS 

We have built several solar models corresponding to different values of Gn, by using 
an up-to-date version of the evolutionary code FRANEC, which includes element diffusion, 
recent opacity tables and modern nuclear reaction rates |TH . 



For a given value of Gn the three input parameters of the code - the mixing length a, 
the initial Helium and metal abundances Y in and Z in are varied until one reproduces the 
observed values of the solar radius [R Q = (6.9598(1 ± 0.01%)10 8 m ], luminosity [L & = 
3.844(1 ± 0.4%)10 26 W] and photospheric composition [Z/X = 0.0245(1 ± 6%)] at the solar 
age [t Q = 4.57(1 ±0.4%) Gy]. 

We remark that the calculated properties of the solar interior are sensitive to the values 
of these observables, particularly to L Q . For this reason, in order to disentangle the effect of 
tiny changes of Gn we require that L Q is fixed to the level of 2 • 10 -5 , i.e. much better than 
the observational accuracy. This precaution is necessary, otherwise the calculated models 
will reflect the changes of L & mixed to the changes of Gn- 

We also remark that astronomical observations fix the product GnM q quite accurately 



G N M Q = (132712438 ±5) 10 12 mV 2 . (1) 

Laboratory measurements of Gat are thus measurements of M & . If Gat is changed M has 
to be varied so that eq. (|TJ) is satisfied. This is an important point for the present discusssion. 

In Figs. |],0,|3] and in table | we present the effect of varying Gn by ±1% with respect to 
the "standard" CODATA-98 value. The following points are to be remarked: 
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i) the changes of pressure and density are of the same order of SGn/Gn, as expected, 
see Fig. [I]. 

ii) On the other hand the induced variation of the (squared isothermal) sound speed 
u = P/p is much smaller. For 5Gn/Gn = 1% one has at most 5u/u = 2 • 1CT 3 , a value 
comparable to the present observational accuracy, see Fig. 0. 

hi) Also the variation of temperature is much suppressed, see Fig. [j| the change of the 
central temperature is about 0.3%, and that of the Boron flux is about 5%, see Table |, in 
both factor three below the present observational uncertainty. 

iv) The depth of the convective zone is altered just by a factor 4-10 -4 and the photospheric 
helium abundance by less than 1%, well below the present observational uncertainties, see 
again Table |I[ 

More generally, we can express the sensitivity of the observable Oi to the change of Gn 
by using scaling parameters 

a log G N 

The calculated fa values are also presented in Table |I[ Changes in Gn at the level of present 
uncertainty (1.5 • 10 -3 ) induce changes of helioseismic observables and neutrino fluxes which 
are negligible in comparison with the respective observational uncertainty, see last row of 
Table | 

Our results do not support the claim of ref. ]I7[ that observational solar data can be used 
to determine the value of Gn with accuracy of few parts in 10~ 4 . Variations of Gn at this 
level induce changes which are much too small in comparison with observational accuracy. 

In summary, present data are sensitive to changes oJGn provided that these are of the 
order of 10 -2 . 



III. INTERPRETATION 

A puzzling situation has emerged. Gn is clearly an important parameter, and its vari- 
ation induces a comparable change on physically relevant quantities such as P and p. On 
the other hand the change of u and T are definitely suppressed. 

The reason for the cancellation is in the constancy of GnMq, as can be easily demon- 
strated. 

Actually the typical scales of solar pressure and density are given by: 

P = G N Ml/R% (3) 
p = M q /Rq (4) 

For a change of Gn which keeps constant the product of GnM q these vary as : 

SP/P = Sp/p = -5G N /G N . (5) 

In fact, this is the pattern shown in Fig. |I[ 

In this approximation, when considering u = P/p the effects of changing Gn cancel, in 
other words u is unaffected by changes of Gn, which explains the much weaker sensitivity 
emerging from the numerical calculations presented in the previous sections. 
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From the perfect gas law, which describes to a good approximation most of the solar core, 
one has P/p = kT/p , so that also T/p is weakly sensitive to changes of Gn- The present 
value of the mean molecular weight p is mainly determined from the initial conditions (Y in ) 
and the solar history and not from Gn- It follows that also T is very weakly sensitive to 
Gn- 

We remark that the constancy o/GnMq is essential for the cancellation. If one computes 
stellar structures with fixed M and different Gn one finds fractional changes of u and T 
that are proportional to SGn/Gn, however these cannot be interpreted as solar models. 



IV. FUTURE PROSPECTS 

One expects that, in the future, the sound speed u will be measured with a higher 
accuracy, possibly to the level of 10~ 4 . Even in this case, however, it will be difficult to get 
significant improvements on Gn, unless one achieves the same level of accuracy in several 
other physical inputs which affect u. For instance, let us consider the effect of variations of 
the nuclear reaction cross sections and of the opacity, k, We remind that these quantities 
are affected by uncertainties of few percent (at least), see [f2T[ -|23| . 



As we have already remarked, the product GnM q is quite accurately fixed, so that if 
Gn increases, the solar mass must decrease according to 5M Q /M Q = —SGn/Gn- A smaller 
Mq, being the radius of the sun fixed by observational data, implies a reduction of number 
densities rii of the various particles present in the sun. 

This implies, a reduction of the nuclear reaction rates (which are proportional to UiUj) 
and, at the same time, an increase of the photon mean free path (which, in the core of 
the sun, is inversely proportional to electron number density). This suggests that suitable 
changes of nuclear energy production rate e and of opacity can mimic changes of Gn- 

Actually, one expects that a combined variation Se/e = —25Gn/Gn and Sk/k, = 
—SGn/Gn should affect the sound speed in the same way as a variation of Gn- This is 
supported by the numerical results shown in Fig. [|. We have built a solar model with 
energy production rate decreased by 2% and with opacity decreased by 1%, and we have 
obtained a sound speed profile which is similar (at the level 0.1% or less) to that obtained 
by increasing Gn by 1%. 

Similar considerations hold for a possible determination of Gat by means of $b measure- 
ments. By using the (3 coefficient shown in Tab. [| one sees that a variation SGn/Gn = 10 -3 
corresponds to change of the 8 B neutrino flux 5&b/&b — 5 • A "global" accuracy (i.e. 
both theorethical and experimental) of the order 5 • 10~ 3 is thus needed in order to distin- 
guish among values of Gn which differ by 10~ 3 . We remind that the present experimental 
and theoretical determinations of the 8 B neutrino flux have uncertainties of the order of 
15-20% per cent.ns hold for a possible determination of Gn by means of <3>b measurements. 
By using the (3 coefficient shown in Tab. |, one sees that a variation 5Gn/Gn = 10~ 3 corre- 
sponds to change of the 8 B neutrino flux 5§b/&b — 5 • 10~ 3 A "global" accuracy (i.e. both 
theorethical and experimental) of the order 5 ■ 10 -3 is thus needed in order to distinguish 
among values of Gn which differ by 10~ 3 . We remind that the present experimental and 
theoretical determinations of the 8 B neutrino flux have uncertainties of the order of 15-20% 
per cent. 
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V. CONCLUDING REMARKS 



We have found that: 



The present uncertainty on the gravitational constant, AGn/Gn = ±1.5 ■ 10 3 has 
significant effects on the profile of density and pressure. 

On the other hand, it has negligible influence on the solar properties which can be 
measured by means of helioseismology and 8 B neutrinos: sound speed profile, central 
temperature, depth and helium content of the convective envelope. 



Our result do not support recent claims [17 that observational solar data can be used 



to determine the value of Gn with an accuracy of few part in 10~ 4 . Present data 
cannot constrain Gn to much better than 10 -2 . Furthermore, even if the sound speed 
measurments will become much more accurate, it will be difficult to get significant 
improvements on Gn due to the approximate degeneracy with variations of other 
physical inputs, e.g. nuclear cross sections and opacity, which are presently known at 
the per cent level or worse. 

ACKNOWLEDGMENTS 

We are grateful to G. Fiorentini for useful comments and suggestions. 



5 



REFERENCES 



[1] S. Degl'Innocenti, W. A. Dziembowski, G. Fiorentini and B. Ricci, Astrop. Phys. 7, 77 
(1997). 

S. Basu, M.H. Pinsonneault and J.N. Bahcall, Ap. J. 529, 1084 (2000). 
W. A. Dziembowski et al. Mon. Not. R. Astron. Soc. 249, 602 (1991). 
J. Christensen-Dalsgaard et al. Ap. J. 378, 413 (1991). 
G. Aardsma et al. (SNO coll.), Phys. Rev. Lett 87, 071301 (2001). 
S. Fukuda et al. (Superkamiokande coll.) Phys. Rev. Lett 86, 5651 (2001) 

F. L. Villante, G. Fiorentini, E. Lisi, Phys. Rev. D 59, 013006 (1999). 

G. L. Fogli, E. Lisi, D. Montanino, A. Palazzo, |hep-ph/ 01 06241 
G. Fiorentini and B. Ricci, Phys. Lett. B. 526, 186 (2002). 
E. Teller, Phys. Rev. 206, 213 (1948) 

M. J. Duff, L. B. Okun and G. Veneziano |physics/0110060| , (2001) : 
P. J. Mohr and B.N. Taylor, Rev. Mod. Phys 72, 351 (2000). 

E. R. Cohen and B.N. Taylor, CODATA Bull. No. 63, pp 1-32. 
J. Luo et al., Phy. Rev. D 59, 042001 (1999). 

J.H. Gundlach and S. M. Merkowitz, Phys. Rev. Lett. 85 2869 (2000). 
T. J. Quinn et al, Phys. Rev. Lett. 87, 111101 (2001). 
I. Lopes and J. Silk, |astro-ph/01h23T0| (2001). 

F. Ciacio, S. Degl'Innocenti and B. Ricci, Astron. Astroph. Suppl. 123, 449 (1997). 
M. Stix, The Sun, Springer- Verlag, Berlin 1991 

J.N. Bahcall, M. Pinsonneault and S. Basu, Astrophys.J. 555, 990 (2001), 
www.sns.ias.edu/ jnb/ 

J.N. Bahcall and M.H. Pinsonneault, Rev. Mod. Phys. 64, 885 (1992). 
NACRE coll., Nucl. Phys. A 656, 3 (1999). 
E. C. Adelberger et al, Rev. Mod. Phys. 70, 1265 (1998). 



6 



TABLES 



TABLE I. SSM predictions, observational errors and sensitivity to Gn- A is the lcr relative 
observational error; 5±\% is the variation (model-SSM) /model obtained when changing Gn by 
±1%; [5 is the scaling coefficient of eq. § ^+0.15% is the variation (model-SSM) /model estimated 
when changing Gjy by +0.15%. 



Observable 


u(O.LR ) 


u(0.6R Q ) 


T c 


<5> B 


Y ph 


Rb 


units 


(Mm/s) 2 


(Mm/s) 2 


10 7 K 


10 6 cm" 2 s- 1 




Rq 


SSM @ 


0.1525 


0.09222 


1.568 


5.15 


0.2437 


0.7140 


A(%) 


0.2* 


0.12* 




18** 


1.4* 


0.2* 


hi% (%) 


+0.18 


-0.036 


+0.31 


+5.7 


-0.67 


+0.04 


(%) 


-0.20 


+0.021 


-0.28 


-4.4 


+0.77 


-0.04 




+0.19 


-0.03 


+0.29 


+5.1 


-0.72 


+0.044 


5+0.15% (%) 


+0.028 


-0.0045 


+0.044 


+0.76 


-0.11 


+0.0065 



*from helioseismic data |l[ 
** from 8 B neutrino data [j|]9] 
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FIGURES 
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FIG. 1. Relative change (model-SSM)/SSM of pressure (upper panel) and density (lower panel) 
as a function of the radial coordinate for a change of Gat by +1% (full line) and -1% (dashed line). 
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FIG. 2. Relative change (model-SSM)/SSM of u = P/p as a function of the radial coordinate, 
for a change of Gat by +1% (full line) and -1% (dotted line). The la (3a) helioseismic uncertainty 
correspond to the dark (light) area. 
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6G N / G N — i 1 0~ 2 
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FIG. 3. Relative change (model-SSM)/SSM of temperature (upper panel) and mean molecular 
weight (lower panel) as a function of the radial coordinate, for a change of Gn by +1% (full line) 
and -1% (dashed line). The vertical bar in the upper panel corresponds to the la uncertainty on 
T c from neutrino measurement. 
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FIG. 4. Relative change (model-SSM)/SSM of u = P/p as a function of the radial coordinate, 
for a change of Gn by +1% (full line) and for a solar model with energy production rate decreased 
by 2% and opacity decreased by 1% (dashed line). The la (3a) helioseismic uncertainty correspond 
to the dark (light) area. 
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